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Abstract. Let Xt be a subordinate Brownian motion, and suppose that the Levy 
measure of the underlying subordinator has a completely monotone density. Under very 
mild conditions, we find integral formulae for the tail distribution P{tx > t) of first 
passage times through a barrier at x > 0, and its derivatives in t. As a corollary, we 
examine the asymptotic behaviour of 'P{tx > t) and its t-derivatives, either as t — > (X) 
or a; — )■ 0+. 



1. Introduction 

The present article complements and extends the results of the recent paper [28], where 
spectral theory for a class of Levy processes killed upon leaving a half-line was developed. 
In a closely related paper ^29j, first passage times were studied for a rather general class 
of one-dimensional Levy processes. In the present article, more detailed properties of 
first passage times are established for processes considered in [2H]: symmetric Levy pro- 
cesses, whose Levy measure has a completely monotone density function on (0, oo). More 
precisely, we prove asymptotic fomulae, regularity, and estimates of the tail distribution 
'P{tx > t) of the first passage time through a barrier at the level x for a Levy process Xf. 

Tx = mi{t>0:Xt>x}, x>0, 

as well as its derivatives in t. Alternatively, the results can be stated in terms of the 
supremum functional Mt = sup^gjo.t] since we have P{tx > t) = P{Mt < x) for all 
t,x>0. 

In [28], a formula was given for generalised eigenfunctions Fx{x) of the transition 
semigroup of the killed process. As an application, the distribution of first passage times 
was expressed in terms of the eigenfunctions Fx{x). The full statement of this result was 
only announced, and a formal proof was given under more restrictive conditions. In the 



present paper, we provide the proof in the general case (Theorem 1.6). The expression 
for the distribution of is then used to find estimates and asymptotic expansion of 
{d / dt)'^P {tx > t). This requires detailed analysis of the eigenfunctions Fx[x). 

The double Laplace transform (in t and x) of P{tx > t) is known for general Levy 
processes since 1957 due to the result of Baxter and Donsker (Theorem 1 in [4j). For 
symmetric Levy processes Xf with Levy-Khintchin exponent "^{C,), 

^-^x-ztpf^^^ > t)dxdt 

1 ^^jATiMillM,,- 



However, the double Laplace transform in (1.1) has been inverted only for few special 



cases. It is a classical result that for the Brownian motion, is the (l/2)-stable subor- 
dinator. An explicit formula for the distribution of was found for the Cauchy process 
(the symmetric 1-stable process) by Darling p.Oj, for a compound Poisson process with 
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^(^) = 1 — cos^ by Baxter and Donsker and for the Poisson process with drift by 
Pyke |32j . A formula for the single Laplace transform for symmetric Levy processes, 
under some mild assumptions, was given recently in ^29j (see Theorem 1.3| below). In the 



development of the fluctuation theory for Levy processes, many new identities involving 
first passage times were derived (see P, [121 EHl E3] for a general account on fluctuation 
theory), including various other characterisations of P{tx > t), at least in the stable 
case, see P[71[TnilIIllIlll5l[l6l[T8l[T9l|26l|271|36]. First passage times and the 
supremum functional Mt play an important role in many areas of applied probability 
([21 13]), mathematical physics ( |20| [25]). and also in potential theory of Levy processes 
(F. "211 [221 [231 [21]). 

The main result of this article is an explicit, applicable for numerical computations 
expression for P{tx > t) for a class of symmetric Levy processes, which includes symmetric 
a-stable processes, relativistic a-stable processes and geometric a-stable processes (in 
particular, the variance gamma process) and many others. More precisely, the following 
assumption is in force throughout the article. 

Assumption 1.1. Any of the following equivalent conditions is satisfied (see Proposi- 
tion 2.13 in [28J): 

(a) Xt is a subordinate Brownian motion, Xt = Bz^, and the Levy measure of the sub- 
ordinator Zt has a completely monotone density. Here Bs is the one-dimensional 
Brownian motion (VarSs = 2s), Zt is a suhordinator (nonnegative Levy process), 
and Bs and Zt are independent processes; 

(b) Xt is a symmetric Levy process, whose Levy measure has a completely monotone 
density on (0, cxd); 

(c) Xt is a Levy process with Levy-Khintchine exponent ^{C) = for some 
complete Bernstein function ipiO- 

We assume that Xt is non-trivial, that is, Xt is not constantly 0. 

Remark 1.2. All explicit formulae and estimates proved in this article are given in terms 
of the complete Bernstein function Translation to the Levy-Khintchine exponent 

^(0 — is immediate, but usually results in less elegant expressions. 

In this article, the term explicit formula is used for an expression involving a finite 
number of (absolutely convergent) integrals, elementary functions and the function ip. 
The complete Bernstein function extends to a holomorphic function on C \ (— oo, 0] (see 
Preliminaries). Sometimes (namely, in the formula for Fx{x)) we also use this holomorphic 
extension of ip. □ 

Our proofs are based on the following two theorems. 

Theorem 1.3 (Corollary 4.2 in |2S])- We have 



f e-«^P(r^ > t)dx = - f 

Jo Jq 




X 

X (1-2) 

pf0^dc]e-'^(''^dX 

foraUt,^>0. □ 

We remark that the above result is proved in [29] for all symmetric Levy processes with 
Levy-Khintchine exponent ^^(0 having strictly positive derivative on (0, oo). 

The transition semigroup pj^^'"^^ (acting on -L^((0, oo)) for any p G [1, oo]) of the process 
Xt killed upon leaving the half-line (0, oo), and its L^{{0, oo)) generator ^(o,oo), are defined 
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formally, for example, in [28j. These notions are only required in the statements of 



Theorems 1.4 and 1.10, and therefore they are not discussed in detail below. 



Theorem 1.4 (Theorem 1.1 in [2H])- For every A > 0, there is a hounded continuous 
function Fx on (0, oo) which is the eigenfunction of Pf^''^\ that is, 



for all t,x > 0. The function Fx is characterised by its Laplace transform: 



X 



exp 



log 



dC 



[1.3) 



for ^ G C such that Re^ > 0. Furthermore, for x > we have 

Fx{x) = sm{Xx + ^x) -Gx{x), (1.4) 

where the phase shift i^x belongs to [0,tt/2), and the correction term Gx{x) is a bounded, 
completely monotone function on (0, oo). More precisely, we have 

''^-nl A^^'°^ *(V)-*(C) 

and Gx is the Laplace transform of a finite measure 7^ on (0, 00). When ip{C,) extends to 
a function ip~^{C,) holomorphic in the upper complex half-plane G C : Im,^ > 0} and 
continuous m G C : Im,^ > 0, and furthermore 7^ 'ip{X) for all ^ > 0, then the 

measure 7^ is absolutely continuous, and 

TT 



X exp 



1 



TT 



1 ..... 

^(A^)-^(e) ^^'^^ 



for ^ > 0. 

We introduce the following two conditions: 



sup 



< 2, 



and, given t^ > 0, 




-to-4'ie) 



dC, < 00. 



:i.6) 



□ 



(1.7) 



:i.^ 



'1 V v^(e) 

Remark 1.5. (a) By Proposition |2.3[ b), for every complete Bernstein function ip the 
supremum in (1.7) is not greater than 2. Condition (1.7), is needed only to assert 



that sup;^>o "^x < 7r/2, and can be replaced by the latter, 

(b) When is unbounded (that is, Xt is not a compound Poisson process) and 

regularly varying of order ^0 at and g^o at 00, then qq, Qoo G [0, 1] and 

.^^oA^\rm/fiO) = 1-00, lim5_oo(el^"(0l/^'(0) = l-Qoo (see ^8]) 



lim, 



Hence (1.7) is automatically satisfied, 
(c) Assumption (|1.8 ) is a rather mild growth condition on ^{$,) for large ^. By 



Proposition 2.3 a), \/W{^^)hK^ ^ Hence, (1.8) is satisfied for all n > 



an all to > whenever '?/'(0 — clog(^) (^ > 0) for some c > 0. When il){C! 



/(log(l + 0), then (|L8|) is equivalent to /r V f'js)/ f{s) e-^''^^''^ds < 00. For 
example, if ■?/'(0 = log(l + log(l + 0), then ([L8| holds if and only if to > 1/2. □ 
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In particular, many processes frequently found in literature, including symmetric a- 
stable processes, relativistic a-stable processes and geometric a-stable processes, sat- 
isfy (1.7) and ( |1.8 ). These and some other examples are discussed in Section [7| 

The following are the main results of the article. The first of them provides a formula 
for P(rj. > t) and its derivatives in t, and it was proved in [28j under more restrictive 
assumptions. The full statement of the theorem (for n = and n = V) was announced 
in |28j as Theorem 1.8. 



Theorem 1.6. // (1.7) and (1.8) hold j or some to > 0, then for alln>0,t> to (t > to 
if n = 0) and x > 0, 



:-ir^p(-.>*) 




(^(A'))"e-*'^(^')FA(a;)rfA. 



In Theorem 4.6 in [29j it is proved that 



:i.9) 



20000 



mm 



' y/tij{l/x^ 



< P(r^ >t) < lOmin 1 



;i.io) 



for all t, X > for a relatively wide class of symmetric Levy processes Xt (similar results 
for many asymmetric processes are also available in [29]). Our second result generalises 
these estimates to derivatives of P{Trc > t) in t, for a restricted class of processes and 



t large enough. For similar bounds with slightly different assumptions, see Lemma 6.2 



Corollary 6.4 Remark |6.5| and Proposition 7.1 



Theorem 1.7. (a) // (1.7) and (1.8) hold for some to > 0, then the distribu- 
tion of Tx is ultimately completely monotone, that is, for each fixed n > 0, 
(-l)"(rf"/c^r)P(r^ >t)>0 fort large enough. 
(b) If g = sup^yoiC\'4^"iO\/'^'iO) < 1 (l-''')y'j then there are positive constants 
Ci(n), C2{n), Csi^n, g) such that 



ci{n) 



< 



' dt^ ^ ' - t-+i/2^/^^(l7^ 



'1.11^ 



t«+l/2^^(l/a;2) 

for n > 0, X > and t > Cs{n, g) / ip{l / x"^) . When C2(n) is replaced by a constant 
C2(n, g), then the upper bound in ( |1.11[ ) holds for all t,x > 0. 

Next, we study the asymptotic behaviour of P{tx > t) as t — i- oo or x — i- 0^. The 
function V{x) is given by an explicit formula; it is the renewal function of the ascending 
ladder-height process (see Preliminaries). The case n = has been previously studied 
in [H]. 



Theorem 1.8. (a) Ifip is unbounded, and (1.7) and (1.8) hold for some to > 0, then 
for all n > and x > 0, 



lim 

t—^oo 



n+1/2 



dt' 



P(r, > t) 



-l)"r(n + 1/2) 



TT 



V{x) 



;i.i2) 



The convergence is locally uniform in x E [0, oo). 
(b) // 'ip is unbounded, regularly varying of order g E [0, 1] at infinity, and ( 1.7[ ) 
and ( 1.8 ) hold for some to > 0, then for all n > and t > to (t > to if n = 0), 

(-l)"r(n + 1/2) 1 



lim 



vW7 



x^ 



dt' 



P(r, > t) 



7rr(l + g) t"+V2 ■ 



;i.i3) 
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The convergence is uniform in t & {to, oo). 

Remark 1.9. The asymptotic behaviour of {d / dt)Y* [t^ > t) as t — i- oo was completely 
described for (possibly asymmetric) stable Levy processes in [13], Theorem 1. In our The- 
1.8[ a), we consider a class of symmetric, but not necessarily stable Levy processes. 



orem 



and derivatives of higher order are included. 

In |13j . Theorem 3, an analogous problem is studied when one-dimensional distributions 
of the Levy process Xt belong to the domain of attraction of a (possibly asymmetric) sta- 



ble law. Again, this partially overlaps with Theorem 1.8 a), but neither result generalizes 
the other one. □ 



Finally, Theorem L4 is complemented by the following completeness result. It was 
proved in |28j under an extra assumption that the operator 11, defined in the statement 
of the theorem, is injective. Using methods developed partially in [29j, we show that 
n is always injective, and therefore the theorem holds in full generality. The present 
statement was announced in 1281 as Theorem 1.3. 



Theorem 1.10. For f E 0^(0, oo)) and A > 0, let 

POO 

n/(A) = / f{x)Fxix)dx. 



(1.14) 
Jo 

Then ^J2/^lY{ extends to a unitary operator on L^((0,oo)), which diagonalises the action 
ofP,^''-\- 



npf '°^V(A) = e- 



-#(A2) 



n/(A) 



/or/eL2((0,oo)). 



Furthermore, f G V{A{o,oo)', L"^) if and only if ip{X'^)Ilf{X) is in L^((0,oo)), and 



nAo,oo)/(A) = -^(A2)n/(A) 



for f e ©(^(0,00);^ 



Remark 1.11. In [28j, the generalised eigenfunction expansion of Theorem 1.10 



the key step in the proof of (a restricted version of) Theorem 1.6 Here, the proofs of 
Theorems 1.6 and 1.10 are independent. □ 



was 



We conclude the introduction with a brief description of the structure of the article. 
In Preliminaries, we recall the notion of complete Bernstein and Stieltjes functions, their 
properties and a Wiener-Hopf type transformation ip ip"^. Some new ideas are developed 
here, e.g. a type of continuity of the mapping ip ip"^ is proved. We also give some 
simple estimates related to Laplace transforms of monotone functions, and introduce the 
renewal function V{x) of the ascending ladder-height process. In Section [3| we prove 
Theorem 1.10 Next two sections contain estimates and properties of 'dx and Fx{x), 



respectively (see Theorem 1.4), which are essential to the derivation of the main results. 



In Section |6| we prove Theorems L6 L8 Finally, some examples are studied in Section [7j 



2. Preliminaries 

2.1. Complete Bernstein functions. A function f{z) is said to be a complete Bern- 
stein function (CBF in short) if 

1 z m{ds) 

Tl Jq+ Z + S S 

where Ci, C2 > 0, and m is a Radon measure on (0, oo) such that j min(s^^, s^'^)m{ds) < 
oo. A function g{z) is said to be a Stieltjes function if 



f{z) =Ci + C2Z + 



(2.1) 



g[zj = — + 02 + - I ^ m{ds), 



z 



TT 



OH 



Z + S 



(2.2) 
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where Ci, C2 > 0, and m is a Radon measure on (0, 00) such that J min(l, s~^)7h{ds) < 00. 

Complete Bernstein and Stieltjes functions are often defined on (0, 00). However, (2.1 ) 
and (2.2) define holomorphic functions on C \ (— oo,0]. We always identify functions on 
(0, 00) with their holomorphic extensions to C \ (— 00, 0]. 

We list some basic properties of complete Bernstein and Stieltjes functions. 

Proposition 2.1 (Proposition 2.18 in [28J, and Corollary 6.3 in [34j). 



(1) Let f be a complete Bernstein function with representation (2.1). Then 



and 



ci = lim f{z), 02 = lim ^ , (2.3) 



m{ds) = lim {Im f{—s + i6)ds), (2-4) 

e— !>0+ 



with the limit understood in the sense of weak convergence of measures. 
(2) Let g be a Stieltjes function with representation (2.2). Then 



ci = lim izg{z)), C2 = lim g{z), (2.5) 

2-S>0+ z^oo 



and 



m{ds) + TTCi6o{ds) = lim {—lmg{~s + ie)ds), (2.6) 

e-!>0+ 

with the limit understood in the sense of weak convergence of measures. □ 

Proposition 2.2 (see [34j). Suppose that f and g are not constantly equal to 0. 

(a) the following conditions are equivalent: f{z) is CBF, z/f{z) is CBF, f{z)/z is 
Stieltjes, l/f{z) is Stieltjes; 

(b) f{z) is CBF if and only if f{z) > for z > 0, and f extends to a holomorphic 
function in C\ (—00, 0] such that Im/(2;) > when Im^ > 0; 

(c) if f , g are CBF and a,c > 0, then also f{z) + g{z), cf{z), f{cz), {z — a)/{f{z) — 
f{a)) (extended continuously at z = a) and f{g{z)) are CBF. □ 

Proposition 2.3 (Proposition 2.21 in [28\). If f is a complete Bernstein function, then 

(a) < zf'{z) < f{z) for z > 0; 

(b) < -zf"{z) < 2f'{z) for z > 0; 

(c) \f{z)\ < (sin(£/2))-i f{\z\) for z e C, \ATgz\ < rr - e, e e (0, vr); 

(d) \zf'{z)\ < (sin(£/2))-i/(kl) for z as m (c); 

(e) 1/(^)1 < c(/,£)(l + \z\) for z as m (c). □ 

Throughout the article, ip usually denotes the complete Bernstein function in the Levy- 
Khintchine exponent of X^. Some preliminary results and definitions, however, are valid 
for more general functions ip. If this is the case, we explicitly state all assumptions on ip 
each time it is mentioned. 

Following [281 129] , we define 

^t(0 = exp^-y^ ^^^dCj (2.7) 

for any positive function ip for which the integral converges. In this case, ip\C} is defined 
at least when Re^ > 0. By a simple substitution, for > 0, 

V.t(0=expf-y^ ,^^2 dCj. (2.8) 



FIRST PASSAGE TIMES FOR SUBORDINATE BROWNIAN MOTIONS 7 

Proposition 2.4 (Proposition 2.1 in [29j). If is a nonnegative function on (0, oo) 
and both and ^/'?/'(^) are increasing on (0, oo), then 

e-'^^/^vW) < i^KO < e'^/'^VW), (2.9) 

where C ~ 0.916 is the Catalan constant. Note that e^^l'^ < 2. 
If, in addition, ip{C,) is regularly varying at oo, then 

lim 4!ffiL = 1. (2.10) 

An analogous statement for ^ — > holds for ipiC,) regularly varying at 0. 

In particular, (2.9) holds for any CBF. Likewise, (2.10) holds for any regularly varying 
CBF. □ 



Tlie estimate (2.9) for CBFs was obtained independently in ^24j, Proposition 3.7, 
while (2.10) for CBFs was derived in [21j, Proposition 2.2. 

Proposition 2.5 (Proposition 3.4 in |28]). Whenever both sides of the following identities 
make sense, we have: 

(a) (l/^)t = (V'^)t = {ij^Y (q e R) and (t/'iV's)^ = V'lV'l; 

(b) {c^ipy = ap^ for c > 0; 

(c) when ?/;(^) = ^, then ^/'t(^) = ^; 

(d) if appropriate limits of exist, then the corresponding limits of ip'^ exist, and 
lim5_o+ ^P\i) = (lim5^o+ lim^-^oo ^-^(0 = (lim^^oo ^(0)'/'- □ 

The first part of the following result was independently proved in [23j, Proposition 2.4. 
Lemma 2.6 (Lemma 3.8 in [28j). //^(O «s a CBF, then also ^t(^) ^s « CBF, and 

= H-e) (2.11) 

for ^eC\R. □ 

Proposition 2.7. Let tpn be a sequence of CBFs. If ipn{0 ~^ "^(0 as n ^ oo for all 

^ > 0, then ipliiO ~^ locally uniformly m ^ G C \ (— oo, 0]. 

Proof. For any CBFs have 



when Re^ > 0. Furthermore, by monotonicity and concavity of ipn, V'n(l) min(l, C^) ^ 
< V'n(l) niax(l, (^^) {( > 0), and therefore there are Ci,C2 > such that 
Cimin(l,C^) < ipniC'^) < C2 max(l, (^^). Hence, by dominated convergence, (2.12) gives 
convergence of V'KO V'HO when Re^ > 0. By Corollary 7.6(b) in [34J, ipniC.) — ?■ V^^O 
for all ^ G C \ (—00,0]. Pointwise convergence of CBFs is automatically locally uni- 
form on C \ (—00,0]: by Proposition 2.3[ d,e), when Arg,^ G (— vr + £:,7r — e), we have 
|(V'D'(OI < C3(£)^i(|el)/lel < C4(£,/)max(|e|-i,l), which proves that are locally 
equicontinuous in C \ (— oo,0]. □ 

As in [28], for a nonvanishing, increasing and different iable function ip with strictly 
positive derivative, we denote 

^^K)= l-'^OMAr «>0,AV«. (2.13) 

This definition is extended continuously by ^/^^(A^) = V^(A^)/(A^'?/''(A^)), and when ip 
extends holomorphically to C \ (— oo,0], also ip\{C,) is defined for ^ G C \ (— oo,0]. For 
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simplicity, we denote = (tpxYl^). By Proposition 2.2[c), if is a CBF, then 

ipx{^) (and hence also is a CBF for any A > 0. 

Although (2.13) is generally used for complete Bernstein functions ip, we will occasion- 
aly need this definition for a negative function ipiO — with g G (—1,0). 

Corollary 2.8. Ifi^iO « CBF, then 

e-''^/-y^^<^l(0<e^'^/'^^^^. □ 

Corollary 2.9. If ip is a CBF, then the function 'ip\{^) is jointly continuous in X > 0, 
^gC\(-oo,0). □ 

The following monotonicity property of plays an important role. 

Proposition 2.10. Suppose that for twice differentiable, nonvanishing and increasing 
functions and ipiC,) we have 



no / -ric) 



< 



Then 



Me) < txie) 



Me) > ±x{e) 



Proof. Integration in ( of both sides of (2.14) yields that 



and by another integration in (, 



C>o. 

< A < e, 

< e < A. 

< ^1 < C, 



(2.14) 



^'(6) - ^'(eo ' ^<^^<^- 

Taking = A^ and ^2 = e> obtain 

Me) _ ^p'{x'){x'-e) ^ ^'ix')ix'-e) _ Me) o<f<x 
M^') M')-M') - i'{x^)-He) M^')' ■ 

The other part is proved in a similar manner. □ 

Finally, recall that tpxi^'^) = '4'{^'^)/{^'^'^'{^^))- This gives the following simple esti- 
mate: 



Me) i^\X''){X^-e) / ^'(A2)max(A2,e') 



< 



A,e>0 



(2.15) 



^a(A2) ^(A2)-^(e) - \M')-M')\ ' 

when ip is a, CBF. 

Let '(/'(O be the complete Bernstein function related with the Levy-Khintchine exponent 
^(^) of the Levy process Xj, \l/(^) = ip{e)- With the notation of this section, for A, ^ > 
we have (see Remark 4.12 in [28], and the proof of Theorem 4.1 in |29| ) 

^ ^^l(Oe-^(^^)c^A; (2.16) 
dC = Argi^litX); (2.17) 



POO r\ PC 

/ e-«^P(r^ > t)dx = - 

Jo Jo 



^9) 



A^-C^ 



log 



A2 + e ^(A2) 

Me 



M^ 
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^m) = ; (2.18) 



,,,, Iv^SMt AV;XA^) 



for the last equality the assumption of the final part of Theorem 1.4 is required. 



2.2. Estimates for the Laplace transform. This short section contains some rather 
standard estimates for the inverse Laplace transform, similar to those used in |29) . 

Proposition 2.11. Let a > 0, c > 1. If f is nonnegative and f{x) < c/(a) max(l, x/a) 
(x > 0), then for any ^ > 0, 

fia) > 



c(l + (aO-ie-<) 



Proof. We have 

pa POO 

e£/(0= / ^e-^''fix)dx+ / ^e-^^fix)dx 

Jo J a 

<cf{a) f ee~^"rfa;+^^ f ^xe-^'^dx 
Jo a Ja 

= c/(a)(l-e-<) + ^(l + aOe-''« 

= c/(a)(l + (aO-V"«). □ 
Proposition 2.12. /// is nonnegative and increasing, then for a, > 0, 

f{a) < e<e^/(0. 

Proof. As before, 

^>C/(0= / ^e-^''fix)dx+ / ee-«V(x)da: 

POO 

> fia) / ^e-^'^dx = f{a)e-^^. □ 

J a 

Proposition 2.13. Let b > a > 0, and suppose that f is nonnegative and increasing on 
(0, b), and f{x) > m for x > b. Then for any ^ > 0, 

1(0') < FI ^7 • 



Proof By Proposition 2A2 applied to /(x)l(o,a)(x) + f{a)l[a,oo)ix) 



f{a) < / e-^'^f{x)dx + / e'^^{f{a) - f{x))dx 

\Jo Ja J 

(POO POO \ 

/ e-^''f{x)dx+ / e-^%fia) -m)dx\ 

= e<e^/(0 + e-(^-")«(/(a)-m). □ 
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2.3. Elements of fluctuation theory. In the sequel we will need the func tion V{x 



which is des cribed by its Laplace transform, CV{^) = Proposition 2.2 a) and 



Lemma 2.6, is a complete Bernstein function. Hence, V" is a Bernstein function, 

i.e. V is nonnegative and V is completely monotone (see By Theorem 4.4 in [29j . 

11 , , 1 

< Vjx) <5 . x>0. 



5 vWM" VWI 



Suppose that extends to a function holomorphic in the upper complex half-plane 
and continuous in the region Im^ > 0. Then by Proposition 4.5 in 



V{x)=bx + - [ Imf-— ^— ^ ^^(l-e-"«)rfe, x>0, 



where b = (limg^o+(^/V'(0))^^^- 

Probabilistically, V{x) is the renewal function for the ascending ladder-height process 
Hg corresponding to Xf. When ip is unbounded, then Xf satisfies the absolute continuity 
condition and V{x) is the (unique up to a multiplicative constant) increasing harmonic 
function for Xt on (0,oo), cf. |35j . 

3. ElGENFUNCTION EXPANSION 
In this section we prove a peculiar integral identity, which (together with the results 



of PH]) yields a short proof of Theorem 1.10 



Lemma 3.1. Let ipiO arbitrary positive, continuously differentiate function on 

(0, oo) with strictly positive derivative, and such that (log(l + '?/'(C^)))/C^ is integrable on 
(0,oo). Then for all ^i, ^2 > ^ , 



^(A2)(A2 + e?)(A2 + ei) 2' 
Proof. As in the proof of Theorem 4.1 in [29j, for ^ > and 2; G C \ (—00, 0] we define 

^ ( ^ r eiog(i+^(aA) 

where log is the principal branch of the logarithm. Clearly, for each ^ > 0, (f{C,, z) is a 
holomorphic function of z, and Lp{E,, z) > when 2; > 0. As it was observed in |29j, when 
Im^; > 0, we have Arg(l + iIj{(^^)/z) G (— tt, 0), so that 

A.,.K..)^-iri^^5(iiipZ£),ce(0../2). 

Jo 6 + 

Hence, for any 6)6 > 0, lm{(f{^i, z)(f{^2,z)) > when Im^; > 0. By Proposition 2.2[ b), 
fiz) = V2(6, z)(p{^2, z) is a CBF of z. 

By monotone convergence, f{z) converges to when z — )■ O"*", and by dominated con- 



vergence, lim2_!.oo f{z) = 1. It follows that the constants Ci, C2 in the representation (2.1 ) 
for the CBF f{z) vanish. Furthermore, for any ,^ > 0, (p{^, z) extends to a continuous 
function ip~^{C,,z) in the region Im^; > 0, given by the formula 

V>-K,^)=exp(^--y_^ <iC 

where log~(z) is the continuous extension of logz to the region Imz < 0. The measure 



m in the representation (2.1 ) for the function f(z) = ^{^1, z)ip{S,2-, z) is therefore given by 
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m{ds) = IT ^ Im(y9+(^i, — s)y9+(^2, — ■s))(is (Proposition 2.l[ a)). By (2.1) and monotone 
convergence, 

1 pim(^+(ei,-s)(^+(6,-s)) 

71 



ds 







-j^ poo 

— lim 



2;Im(v9+(^i, -g)v?+( 6, -s)) ds 

s 



(3.2) 



z + s 

lim {(p{^i, z)(p{^2, z)) = 1. 



Note that if ipiO bounded on (0, oo) and s > sup^^QipiOy then (p~^{^,—s) is real, so 
that the integrand in (3.2) vanishes for such s. Hence, we may substitute ipiX"^) for s 
in (3.2) to obtain that 



A^'(A') Im 



^+(6,-^(A^))v^+(6,-^(A^)) 



dX = 1. 



Finally, in the proof of Theorem 4.1 in [29j (formula (4.3) therein) it is proved that 
Since Im((A + ^ii)(A + ^2^)) = A(^i + ^2), the lemma is proved. 



Remark 3.2. Lemma 3.1 was conjectured in a preliminary version of |28] 



□ 
□ 



Proof of Theorem \l.lC\ In [28J, Theorem 1.10 was proved under an extra assumption that 
the operator 11 defined in (1.14) is injective. Below we show that this condition is always 
satisfied. 



Let eg(x) = e-«^l(o,oo)(a;) (^ > 0). By ( |2l8| , for A > 

/■oo 

ne^(A) = / e-^''Fx{x)dx = CFxiO = 
Jo 



7A(A2) A2+e^^^' 



By Lemma 3.1 



- AV(A2)V;1(6)^1(6) 
((0.-)) J, ^(A2)(A2 + ef)(A2 + e|) 



dX 



The functions (^ > 0) form a linearly dense set in L^{{0,oo)). By approximation, 
y^/yr n is a unitary operator on i^^((0, cso)). In particular, 11 is injective. □ 

4. ESTIMATES OF 'dx 

As it will become clear in the next section, upper bounds for the phase shift i^x (which 



is defined by (1.5)) are crucial for applications of Theorems 1.4 and 1.10 when -^x is 
close to 7r/2, estimates of the eigenfunctions Fx{x) are problematic. Recall that by ( |1.5[ ) 
and (|2l7l). 



^x = ArgipliiX) 



Jo 1 - V'a(A2) 



A > 0. 



(4.1) 



Here is a complete Bernstein function in the Levy-Khintchine exponent of Xt. By 
a substitution ^ = 1/s in the integral over (1, 00), one obtains that (see Proposition 4.16 
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in [28]) 



1 



1 



log 



A > 0. 



(4.2) 



When "ipiO — ^"^^ (0 < a < 2), that is, when Xt is the symmetric a-stable process, we 
have = (2 — 0)11/8 (Example 6.1 in [28]). In the general case, < -^x < 7r/2, and 
''^A < (^/2) sup^>o(^^a(0/'^a(0) (Proposition 4.18 in [28J). This estimate is insufficient 
for our needs, and it is significantly improved below. By (4.1) and Proposition 2.10, we 
have the following monotonicity result 



Note that ip need not be a complete Bernstein 



function in Propositions 4.1 and 4.2 



Proposition 4.1. Suppose that ip and ip are arbitrary twice- differentiate , increasing 
and nonvanishing functions (not necessarily complete Bernstein functions). Define -t^x 

We assume that the integrals in the 

-rio/fio < -rio/fio for 



by (4.2), and define -t^x in a similar way, using ip. 



definitions of-dx and -dx o-re convergent and finite. If 

all ^ > 0, then i9x < ^\ for all A > 0. □ 

Taking ^'(0 = (corresponding to the symmetric a-stable process) and using the 
fact that -i^A = (2 — a)7r/8, we obtain a global estimate for -t^x- It turns out that although 
this method works for a G (0,2], it generalises easily to a G [—2,2]. First, we recall the 
following result, which was given in Example 6.1 in |2H] for g G (0, 1]; the proof, however, 
works for any g ^ 0. 

Proposition 4.2 (see Example 6.1 in [28]). Let g ^ and^{i) =^s if g>{], ^{^) = 

if g<Q. Define §x by (|42|). Then ^9;, = (1 - ^)7r/4. □ 

From now on, ip is again the complete Bernstein function such that ^'(0 = "^(^^ 
the Levy-Khintchine exponent of the process Xj. 

Proposition 4.3. We have 



IS 



inf 



TX 



4 



sup 



TT 

4' 



(4.3) 



Proof. Denote the supremum in ( |43| by 1 - ^. By Proposition Qb), g G [-1,1]. 
Suppose that g is non-zero, and let ip{C) = if > 0, ip{C) = — if g < 0. Observe 
that ip is increasing and —^ip" {0 / iC) = (1 ~ f?); so that 



C>o. 



^'(0 - v^'(C) 

The upper bound in ( |4.3 ) follows by Propositions 4.1 and 4.2 When ^) = 0, we simply 
consider ip{C) = —C,^^ and let e — )■ O"*". Finally, the lower bound is proved in a similar 
manner. □ 

We conclude this section with a local estimate of 'dx, which depends solely on ip{X'^), 
ip'{}?) and %p"{\^). This result is used only in the construction of an irregular example 
in Subsection |7.6[ We need the following technical result. 



Proposition 4.4. For a > 0, 

1 



TT 



, 1 + a^x^ , 
log -— — ^ dx 
1 + 



arctana. 



For he (0,1), 



vr 



1 -x2 



log(l - 6^(1 _ x^))dx 



(arcsinJ))' 



TT 
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Proof. The function log((l + a^x^)/(l + a^)) is holomorphic in the upper half-plane with a 
branch cut along {ia~^,ioo). By an appropriate contour integration and limit procedure 
(we omit the details), 



, 1 + a^x^ , 
log — — ax 



dx 



-27r 



x^ 



l/a 



y 



, 1 + a^x^ , 
log — — — dx 
I + 

dy = — 27rarctana. 



x^ 



For the second equality, we have by Taylor expansion and beta integral. 



poo T ^ i2n pi 

/ \og{l-b\l-x'))dx = J2— 

oo ,2n /•! °° 
n=l "^'J n=l 



X 



2\n-l 



dx 



T{n)T{l/2)b' 



2n 



2nr(n + 1/2) 



(arcsinfo)" 



see e.g. p. 108 in |3T] for the last identity. 
Proposition 4.5. We have 



□ 



TT 



■(^A < arcsin 



and 



> arcsin 



'A2|V'"(A2)| 



2V''(A2 



+ arcsin 



A2^/>(A2)|V'"(A2)| 



( 



arcsm 



A4|^/>"(A2)| 



2V''(A2)(^(A2) - AV(A2)) 

2 



2(^(A2) - AV(A2)) 



It is easy to see that these bounds are sharp for ipiO — Ci^/(c2 + C,) and ipiO — 
(and, in fact, these are all CBFs ip for which equalities hold), see examples in Section [7| 

Proof. Note that ^pxiO) = 1 and V'a(A^) = ^p{X^)/{X^^'{X^)). Let = ^pxiX^) - 1- Since 
ipx is concave, we have ipxiX'^C) > 1 + a^C^ for C < 1; and ^x{X'^0 < 1 + a^C^ foi' C > 1- 



Hence, by Proposition 4.4 

< 



1 ^a(A^C^) 



1-e""" V^a(A2) 

^logi±4!.C^arctana. 
^ 1 + 



1-e 



Since arctana = it/2 — arcsin(l/ (1 + a^)^/^), this proves the upper bound. 

For the lower bound, let aj = ^/'a(A^) = 1 + and al = A^^/'^(A^)/-?/'a(A^). By a short 
calculation, a"^ = A^|'?/'"(A^)|/(2?/''(A^)). Since ^/'a(0 is concave, its graph lies below the 
tangent line at ^ = A^, so that 

^a(a^C')<^a(a^)-aV;(a^)(i-C') 

= alii- 4(1 -e)). 
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A similar argument for the CBF ^/{ipxiO ^ 1) ^ short calculation give 



^2 „2 „2„2 



Qj ^1^2//- \2^ 



which for ^ = A^/C^ (C > 0) reduces to 

V^a(AVC')>1 + 



4a-2 



a^ + ajalil-C) ' 



Therefore, by (4.2) 



TX 



> 



1 



1-C 
1 



1 , M^va 

2 log dC 



log 



l-a-'al{l-C 



□ 



ttJ, 1-e ^ {1 - alii - - a-^alal{l - e)) 

Triple application of Proposition |4.4| yields 

TTi^A > (arcsina2)^ + (arcsin(aia2/a))^ — (arcsin(a2/a))^. 

5. Properties of Fx{x) 

Some basic estimates of the eigenfunctions Fx{x) have already been established in 
However, for the proof of Theorems 1.6 -1.8, more detailed properties of Fx{x) are re 
quired. Recall that by (1.3) and (2.18), 



A 



and by (1.4) 



Fxix) = sm{Xx + ^x) -Gxix), 



A > 0, Re^ > 0, 



A > 0, X > 0. 



(5.1] 



The phase shift i)x was studied in detail in the previous section. Recall that the correction 
term Gx{x) is a completely monotone function, Gx{x) = £7a(x), and 7^ is a finite measure 
on (0,00). In many important cases, 7^ is given explicitly by ( |1.6 ) and (2.19); however, 
in this section these identities are not used. We extend the definition of Fx{x) by letting 
Fx{x) = for X < 0. 

It was proved in Proposition 4.22 in [28j that 



m(0 < 



A2 + e2 - A + r 



< 



A,e>0. 



(5.2) 



By combining (5.1) and Proposition 2.4 we obtain the following more detailed estimate. 



Corollary 5.1. We have 

X 



'Me) 



A2+e v^a(a2) 



< CFxiO < e 



2C/n 



X 



'Me 



A2+e v^a(a2) 



/orA,e>0. 



(5.3) 

□ 
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Holder continuity of F\ was already studied in Lemma 4.24 in [28]. However, when 
ip does not have a power-type growth at infinity, Fx fails to be Holder continuous. On 
the other hand, uniform estimate for Fx{x) for small x > is crucial in the proof of 
Theorem 1.6 We generalize the estimates of |2B] in the two lemmas below. 

Since Gx{x) is completely monotone, we have G\{x) < and G'l{x) > for A, x > 0. It 
follows that the function Fx{x) = sm{Xx+^x)—Gx{x) is increasing on [0, {71/2— {}x)/ X] and 
concave on [0, (tt — ■})x)/X] (this explains the role of upper bounds for i)x)- Furthermore, 
< Gx{x) < sm^x (Lemma 4.21 in [28j), and so -1 - sin^^A < Fx{x) < 1. 

Lemma 5.2. When A > and < Xx < (7r/2 — ^x)/'^, we have 



(5.4) 



n/2-^xJ V ^a(A2 



Proof. Let A>0, 6=(7r/2 — i)x) / X and < a < 6. By concavity of the sine function and 
convexity of Gx, for any x > a we have 

sin(Ax + -(^a) ~ sin-i^A ^ sin(Aa + {}x) — sint^A 
X ~ a ' 

sint^A - G'a(x) ^ sin^A - Gx{a) 
X ~ a 

Hence Fx{x) < {x/a)Fx{a) when x > a. This and monotonicity of Fa on [0,6] yield 

Fx{x) < Fx{a)max{l,x/a) 



for all X > 0. By Proposition 2.11 for ^ > we have 

^CFxiO 



Fx a ) > 



1 + (a^-^e-^' 
For ^ = 1/a, this gives 

e + 1 a 

On the other hand, -Fa(x) is increasing in x G (0, b), and clearly -Fa(x) > —2 for all x > b. 
Proposition 2.13 gives 

e^^^CFxiO + 2e-(^/2-^A-Aa)5/A 



Fa a < 



1 _ ^-{n/2-^x-Xa)i/X 



for all ^ > 0. Let C, = 1/a. Using the inequality 1 — e ^ > s/(s + 1) (s > 0) in the 
denominator and e~^/^ < s (s > 0) in the numerator, we obtain that 

Fxia) < . . . ( ^^^MIM + 2e-(-/2-'^.-Aa)/(Aa) 



1 _ (,-{-n/2--dx-\a)/{Xa) 



a 



e f CFx{l/a) _^ 2e-(-/2-^A)/(Aa) 



1 _ g-(7r/2-i?A-Aa)/(Aa) 

^ e(7r/2 - ^a) f CFxil/a) 2Xa 
- 7r/2 - ?9a - Aa V a n /2 - ^ 
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The above bounds for F\{a) combined with (5.3) give 



e + i P + l/a^V V'aIA^) 



Assume now that a E (0,6/2]. Then a < b/2 = {tt/2 - ^9a)/(2A) < 1/A, and therefore 
< + l/a2 < It follows that 



2(e + l) V ^a(A2) 



e(7r/2-^/ / / ^,(l/a2) 2Aa 
- n/2 -^,-\a V Va(A2) + ^/2 - 

Since V'aIO is increasing on (0, oo) and A < 1/a, we have 1 < >yilJx{l/a'^)/ilj\{X^). Finally, 
Aa < A6/2 = (7r/2 - so that (7r/2 - ^x)/{t^/^ - ^\ - Aa) < 2. Formula Q 

follows. □ 

Since Fx{x) = sin(Ax + -t^x) — Gx{x) for a completely monotone Gx, the modulus of 
continuity of Fx is described by the behavior of Fx{x) for small x. More precisely, we 
have the following result. 

Lemma 5.3. We have 



,r..N r. / M • /30A|x-i/| ipxil/lx-yl'^) A 

\F.(.) - F,iy)\ < mm [-^J^ f ' ■ (5.5) 

for all X > and x,y > 0. 

Proof. The inequality \Fx{x) — Fx{y)\ < 4 is clear. For x,y > and A > 0, 
\Fxix) - Fxiy)\ < \Gxix) - Gx{y)\ + X\x - y\ 

< \Gxi\x-y\)-GxiO^)\ + X\x-y\ (5.6) 

< \Fx{\x-y\)\ + 2X\x-y\. 

The same inequality is obviously true also when x = or ?/ = 0, so from now on we 
assume that x,y >0. 

Suppose that X\x - y\ < {n/2 - dx)/2. Then, by and (Q, 



As in the proof of Lemma 5.2, 1 < y^^pxi^/l^ ~ 2/P)/V^a(A2) (^indeed A < l/\x — y\, and 
is increasing), so that 



4e , ^\ , hpxil/\x-y\^) 



\Fx{x) - Fxiy)\ < ^26^+^'^/'^ + ^^^^^ + 2^ A|x - y| 



Finally, the parenthesised expression does not exceed (vre^+^c/Tr _j_ 4g _|_ 7]-)/(7r/2 — i^a) 
and (5.5) follows. 
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It remains to consider X\x — y\ > (7r/2 — 'd\)/2. We claim that in this case, the 

mmimum m d^Sj ) is equal to 4. By Proposition Qa), /(O = ^/^a(0 

is a CBF. Hence, 



f{x)/f{a) > min(l,x/a). We obtain 




AV^a(A2) 



2^ 



V {l/\x-y\^)/M^/\x-y\ 
min(A|x — y|, 1) 1 

This proves our claim, and the proof is complete. □ 

The behavior of Fx{x) when x 0^ (with fixed A > 0) was studied in Lemma 4.27 
in [2S]: if ip is an unbounded CBF regularly varying of order a G [0, 1] at infinity, then 
Fx{x) is regularly varying of order a at 0, and 

hm (^/W/^)Fx{x)) = A > 0. (5.7) 

The next result shows that as A — ?■ 0+ with fixed x > 0, Fx{x) behaves as Xy^'ip'(X'^)V{x), 
where V{x) is the renewal function of the ascending ladder-height process, CV{^) = 
^/^^(O (see Preliminaries). 

Proposition 5.4. Suppose that i/j is an unbounded CBF. The function Fx{x) is jointly 
continuous in X > 0, x >0. Furthermore, z/ lim supA_^o+ ''^a < 7r/2, then 

F^(x) 

lim = V(x), x>0, 

^^0+ Xy^WW) 



and the convergence is locally uniform in x > 0. In other words, Fx{x) / (X^/^tf/JX?)) 
extends to a continuous function in X,x > 0. 

Proof. Roughly speaking, we use estimates of CFx to show continuity of e~^Fx(x) in 
A > 0, and then estimates of Fx to replace convergence by locally uniform convergence. 



Let A„ — )■ A > 0. By Corollary 2.9, £Fx„(l + i^) converges to CFx{l + i^) pointwise 
for e R. By (5.2) and dominated convergence, CFx„{l + i^) converges to CFx{l + i^) 
in L^(R). Note that CFx{l + iC,) is the Fourier transform of e~^Fx{x). By Plancherel's 



theorem, e~^Fx„{x) converges to e~^Fx{x) in L^(R). By (5.5), the sequence Fx„{x) is 
equicontinuous in x (here we use the assumption that ip is unbounded), and hence it 
converges to Fx{x) locally uniformly in x > 0. The first part of the proposition is proved. 
As A —7- 0^, the functions (A^/'0(A^))'?/'a(O converge pointwise to the CBF C,/ip{0- 



Therefore, by Proposition 2.7, ^JX^ /ipiX"^) V^a(0 converges to i/ip'^{C) (^ G C \ (— oo,0]). 
We conclude that £Fa(1 + ii)/{X^^'{X'^)) converges to CV{1 + i^) (^ G R). Again, 
we obtain i^^(R) convergence of e^^Fx{x) / {X\J ip' {X"^)) to e~^V{x). The proof will be 
complete if we show that the family of functions Fx{x) / {X^J1l)'{X'^)) is equicontinuous as 
A^0+. 



By (5.5 ), we have 

\Fx{x) - Fx{y)\ 

x,/WW) 



< 



30 X — y\ 


li)x{l/\x-i 


1?) 


7r/2 - ^A V 
30 


V^'(A2)V;a(A2) 


/ 1-A2|. 


X — 


vr/2 - ?9a V 


i){l/\x - y 


2)-^(A2) 



(5. 
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Fix Ao > small enough. By the assumption, tt/2 — 'd\ is bounded below by a pos- 
itive constant for A G (0, Aq). It follows that for A G (0, Aq) and x,y > satisfying 
\x — y\ < 1/Ao, the expression on the right hand side of (5.8) is bounded above by 
c/ ^yip{l/\x — yl)"^ — '0(Ao) (with c depending only on ip and Aq). This upper bound does 
not depend on A G (0, Aq), and since is unbounded, it converges to as |x — 1/| — )■ 0. □ 



6. SUPREMA AND FIRST PASSAGE TIMES FOR COMPLETE BERNSTEIN FUNCTIONS 



Below we prove Theorems 1.6 -1.8, We remark that the condition (1.7) is used only to 
assert that sup^^^g'^^A < 7r/2, and can be replaced by the latter condition. 



Proof of Theorem L6. First we consider n = 0, that is, we will show that 

2 



P(r, > t) 




(6.1) 



We claim that the integral in (6.1) converges, and that e~^^ e"'^^^^'^'' F\{x) \J ip'^X^) /ip{X^) 
is jointly integrable in A, x > for any > 0. Assuming this is true, the proof is quite 
straightforward. Indeed, by Fubini, the Laplace transform (in x) of the right hand side 
of (6.1 ) is then 




By Theorem 1.3 (see (2.16)) and Theorem 1.4 (see (2.18)), this is equal to the Laplace 
transform (in x) of P{Tx > t), and the result follows by the uniqueness of the Laplace 
transform. Hence it is enough to prove our claim. 

Let t > tQ, Q = sup^^^Q-i^A and Aq = (7r/2 — Q)/{2x). By the assumption (1.7) and 
Proposition 4.3, O < tt/2, and hence Aq > 0. Note that XqX < n/A < 1. Since |Fx(x)| < 2 
and ^^/''(A2)/V;(A2) < 1/A, we have 




< 2 



-#(A2) 



A 



dX + 2 



'min(Ao,l) 

< 2max(0,-logAo) + 2 




^(A2 



which is finite by assumption (1.8). We now consider A G (0,Ao). By Lemma 5.2 
have 



we 



< 



,/^e--^(^-)FA(.)<cUe)A.j ^^^''^^^^ iM 
V ^(A2) J^x[x)<c,[U)Ax^ ^(A2)^.(A2) 



Furthermore, by (2.15) (recall that A < Aq < 1/x) 



Ax 



ltlj'{X^)Ml/x^) 



< 



A^'(A2 



V'(A2)^a(A2) - ^V'(A2)(^(l/x2)-7A(A2)) 
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Integration and substitution z = ipi^"^) give 







- " 'Jo Vv^(A2)(V'(l/x2)-^(A2)) 

ci(e) /•'^(^o) 1 rf^<£i(^ 



2 Jo ^zm/x^)-z) 
We conclude that 



^2^e-*^("')|F,(x)|rfA < C2(V^,t) + 2max(0,logx), 



which shows that the integral in (6.1) is absolutely convergent. This also shows that 



^0 



as desired. The proof in the case = is complete. 

For general n > 0, one uses dominated convergence to prove that the derivative can by 
taken under the integral sign; we omit the details. □ 

We now find upper and lower bounds for the derivatives in t of P{tx > t). Note that 
the estimates of P(r2. > t) are covered by [22] for a wider class of Levy processes. We 
begin with a simple technical result. 

Proposition 6.1. (a) Let 7(0; x) = e~'^s"'~^ds (a,x > 0) be the lower incomplete 
gamma function. Then 

min(l,x") min(r(a + l),a;") 

<7(a;x)< , a,x>0. (6.2) 

ae a 

(b) We have 

"^e-'s""^, min(2r(a) + 7ra"e-",2(l/a + , , 
ds < ^ — — ^ ' —, a,x>0. (6.3) 

VX - S y/2x 

Proof. For the part (a), we simply have 

, ^ r -s a-1 , ^ 1 r'"^''"^ min(l, x'^) 

7(a;x) = e s ds > - s ds - 

Jo e Jo 



e In ae ' 



and 

7(0; x) = I e'^s^'^ds < min ( r(a), / s'^'^ds ] = min ( T{a 



a 



JO \ JO 

To prove (b), we split the integral into two parts. First, 

fx/2 /~2' rx/2 



X — ~S V 



ds< \l- I e-'s^-^ds 







2 2 f x"- 

7(0; x/2) < \ — min I T(a), — 
X \ X \ a 
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Next, s'^e"^ attains its maximum at s = a. Hence, 

fX s a-l fX 1 

/ %== ds < a^e-"^ / , 

Jx/2VX-S A/2VS2(X-S) 



ds 



x/2 ^Js{x - S) 



ds 



Finally, 



It follows that 



x/2 



y/x - S 



ds < X 



a-l 



x/2 



\/x - S 



ds = v^.t'^-^/^ 



X g-.^a-l 







,2 min(r(a + l),x'') 
ds < \ — - + min 



X — s 



< 



2 /min(r(a + 
X \ a 



+ min 



, X 



which gives (6.3). 



□ 



Lemma 6.2. Suppose that (1.7) and (1.8) hold for some to > 0. For n > 0, t > to 

(t >to if n = 0) and Aq > 0, denote (of (1.8)) 



4 roo 



Ao 



Let 6 = supg^Q^^A md Ao(x) = (7r/2 — 0)/(2x). If t > to (t > to if n = 0) and x > 0, 
then 



2\\n 



Ci{n,Q) min {^Pil/x')) 



1 



t"+l/2^^(l/a;2) 



+ In{t, Ao(x)) 



(6,4) 



2\\n 



< C2(n,e)min {^Pil/x')) 



t"+l/2^^(l/a;2 



In(t, Ao(x)). 



Remark 6.3. (1) When (1.8) holds for some to > 0, then clearly /o(to,Ao) < c>o for 
any Aq > 0. Furthermore, for any n > 0, t > to and Aq > 0, we have e~(*~*o)'^(''^ ) < 
c(n, Ao)(V'(A^))"' for A > Aq. Hence, (1.8) implies also that /„(t,Ao) < oo for any 
n > 0, t > to and Aq > 0. 
(2) When t > to, then /„(t,Ao) < e-(*-*«)^(^o)/„(to, Aq). Hence, In{t,\o), if finite, 
decays exponentially fast as t — ?■ oo. Therefore, /„,(t, Aq) can be considered as the 

□ 



'error term'; see Theorem 1.8 



Proof of Lemma \ 6. ^ By the assumption ( |1.7[ ) and Proposition 4.3, O < 7r/2. Let k = 
(7r/2 — 6)/2 G (0,7r/4), so that Aq = Ao(a;) = k/x. Denote the integrand in (1.9) by 

/n,t,x(A), 



fn,t,xW 



1^ (^(A^))"e-^(^^)F.(a:) 
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By Lemma 5.2, for A G (0, Aq) we have 

fn, 



,l-2CAx„,,//x2 



2(e + l) 
=i-2CAA7//rA2l 



v/^(l/x2)-^(A2) 



" 2(e + l) 

Since Xx < k and '?/'(Aq) = ip{k'^/x'^) > k'^ip{\/x'^) {ip is nonnegative and concave, and 
k < 1), we obtain (with z = ■?/'(A^)) 



Ao 



fn,t,xWd'>^ > 



> 



4(e + l)vWW) -/o 
4(e + 



7(n + l/2;fc2tV^(l/x2)), 



where 7 is the lower incomplete gamma function. By (6.2) 



mm 



in(l, (A;2tV(l/a;2))"+V2) 



k 



4(e + l)t«+i/2^^(l/x2 
In a similar way, 

/n,,.(A) < (^26^+2^/- + 

and so 



(n + l/2)e 



VV^(1/X2)-V^(A2) 



(6.5) 



^/'{l/l^^) ^n-1/2 -i2 



\/t?(l7^ 



(is. 



By (|6]3|), with ci{n) = 2r{n + 1/2) +7f{{n + l/2)/e)"+i/2 ^nd C2{n) = 2{l/{n + 1/2) + 1), 

r-Ao 

(6.6) 



e\ mm 



in(ci(n), C2(n)(tV^(l/x2))"+i/2) 



t"v^2tV(l/a:2 



Hence, we found a two-sided estimate for the integral of fn,t,x{^) over (0, Aq). The integral 
over (Ao, 00) is highly oscillatory, and therefore difficult to estimate. For this reason, we 
are satisfied with a simple bound obtained using the inequality |Fa(s)| < 2, 



Ao 



\fn,t,xW\dX < - / mX')r--'/'./i7i)^)e-'^^''UX 

^ -'Ao 



(6.7) 



In(t, Ao). 



The lower bound in ( 6.4[) is a consequence of (6.5) and (6.7), and the upper bound in (6.4) 
follows from (6.6) and (6.7). □ 
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We remark that in the statement of the lemma, we can take 



Ci(n,e) 



- 7rVl6 (7r/4 - 0/2)2"+i ^ {■K/2-Qf"+^ 



27re(e + l)(n + 1/2) 



17(n + l/2)22"+3 



3^2 



TT 



7r/2 - 



{2V{n + 1/2) + 7r((n + l/2)/e)"+i/2) < 



15 ri! 
7r/2 - e ■ 



Note that Ci(n, O) decreases with G, while C2{n,Q) increases with B. The notation of 
Lemma 6.2, namely In{to,)\o), Ci{n,Q) and C2(n, G), is kept in the remaining part of 
the section. 



Corollary 6.4. Let e>0,n>l,tQ>0, Xo>0. If the conditions (1.7) and (1.8) hold 
true, then there are positive constants Cs^n, 0), C4(n, 9), C^^e, n, 6, In(to, Ao(a;o))) (here 
6 = supg^o^^A and \o{x) = (7r/2 — 9)/(2x)j such that 



^'»'^' < (-1)" P(r. > *) < ^'<"-Q' 



(6.8) 



w/ien X G (0, Xq] and 



t > max (1 + e)tQ 



1 C5(£,n, 9,4(to, Ao(Xo))) \ 

V'(l/X2)' (V'(l/x2))l+^ ' 



Proof. This is a combination of Lemma 6.2| and Remark 6.3[ As in the proof of Lemma 6^ 



we let k = (vr/2 — 0)/2, so that Aq = Ao(a;) = k/x. Recall that since A; < 1, we have 
?/'(fc2/a;2) > k'^ip^l/x'^). It follows that for n > 1, £ > 0, a; > and t > {1 + e)tQ, we have 



In(t, Ac 
Inito, Ao) 



< e 



-fc2£/{l+£)tV)(l/x2) 



< 



02(6, n,Q) 



(tV^(l/x2))"+l/2+n/e • 



Fix A>0. When )f:(^(l/a;2))i+^ > we obtain 

€2(6, n,Q) 



/n(t,Ao) < 



A'^t"+V2^^(l/a;2; 



-^n(^05 Aq). 



Hence, by (|6.4 ), if ttpll/x'^) > 1, we have (with the constants Ci{n,Q) and C2(ti, O) of 
C2(e,n, 0)/„(to, Aq) 



Lemma 6.2 



Ci(n,e) 



1 (i" 
' < (-1)" i- P(t, > t) 

t"+V2^^(l/a;2) - ^ ^ rft" ^ " 



< C2(n,e) + 



C2(£,n, 9)/n(to, Aq) 



t«+l/2^^(l/a;2 



When X G (0,xo], then /„(to, Ao(a;)) < /„(to, Ao(xo)). Hence, (6.8) holds if A > 
{C^{e, n, 6, /„(to, \o{xq)))Y/^ for some C5. □ 



Proof of Theorem ] 1.11 Part (a) follows directly from Corollary 6.4 For part (b), suppose 
that \'ip"{0\/'^'{0 < ^/^ foi' some g G [0, 1) and all ^ > 0. By Proposition [4!3| with the 
notation of Lemma [6. 2| we have G < gn/A. 

Integrating the inequality —'ip" (0 / i'' iO — qUi '^^ obtain 'ip' / "ip' {^2) < (^2/^1)^ if 
< ^1 < ^2- Integrating this again in ^1 gives ip{i) /ip' {C) < ^/(l — g). Hence, for all 
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In{t,Xo) < 



-#(A2) 



Ao 



(^(A'))"-^/'rfA 



2 



-*^(^')(t^(A'))"-^(2AtV^'(A'))rfA 



Ao 



r(n;t^(A^)), 



where r(a; 2;) is the upper incomplete Gamma function. In particular, In{t, Aq) is finite, 
and (1.8) holds true. 

{n/2 - e)/2 G (7r/8,7r/4), and take Aq 



Let k = (7r/2 - e)/2 G (7r/8,7r/4), and take Aq = Ao(x) 
Recall that since k < 1, we have iP{Xq) > k'^ip{l/x^). Hence, 

2 T{n;kHij{l/x^)) 



k/x, as in Lemma 6.2 



/n(t,Ao) < 



2 yt^<l7x2)r(n; (7r/8)2t^/>(l/x2)) 



(6.9) 



The constant C3(n, f?) in the statement of the theorem is so chosen that for all s > C3(n, g), 



vrVl - g 



min(Ci(n, 7r/4), C2{n, 7r/4)) 



with the constants Ci{n, Q) and C2(n, G) from Lemma 6.2, Then, by Lemma 6.2, when 
tipll/x"^) > C3(n, g) we have 

2t"+V2^^(l/a;2) - ^ ^ dt" ^ " ^ - 2t"+l/2^^(l/^2)' 



where 
6.2 prove that 



as desired. Finally, by (6.9), we have /„(t, Aq) < c{n, g)/{t^^'^/'^ 
c{n,g) = 2n^^(l — f))"^/^ sup^>Q(s^/^r(n; (77/8)^5)). This and Lemma 
the upper bound in (1.11 ) holds for all t,x > 0, but with a constant depending on g. □ 

Remark 6.5. (a) The strict inequality ^ < 1 is essential to the proof. The case g = 1 



is much more complicated; see the example in Subsection 7.3 



(b) Theorem 1.7[ b) can be easily generalised to the case when |V^"(OI/V''(0 — 0/^ 
only for ^ G (0,e) U (1/e, 00), with e G (0,1) fixed. In this case the constant 
C3 = 03(11, g,e) depends also on e. We omit the details. 

(c) The upper bound in (1.11) is certainly not optimal when tipll/x'^) is small. This 
is due to essential cancellations in (1.9). □ 



Proof of Theorem 1.8(a). We use the notation of Lemma 6.2, and take Aq = Ao(a;) 
(7r/2 — 0)/(2x) for a fixed x > 0. Let fn,t,x{X) be the integrand in (1.9). We have 



Ao 



/ /n,t,x(A)rfA 



Ao 



Fx(x) 



n-l/2^-tV(A2) 



By Proposition |5.4[ limx^o Fx{x)/{X^/WU^) = V{x), and therefore Fx{x)/{X^/ip'{X'^)) 
extends to a continuous function of A G [0, Aq]. Let 



fitidX) = e+'/\^{x')) 



n-l/2„-t^(A2) 



AV^'(A2)l[o,Ao](A)rfA. 



As t — )■ 00, the density function of fit converges uniformly to on [e, Aq] for every e > 0. 
Hence, fit{dX) converges weakly to a point-mass at 0. Furthermore, by a substitution 
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\M 



7(n + l/2;tV^(A2)) 



and hence \\^t\\ converges to r(?2 + l/2)/2 as t — oo. It follows that 



lim I 



Ao 



r(n + i/2) 



V{x) 



Finally, by (6.7) and Remark 6.3 

4 1-°° 



lim 



n+1/2 



fn,t,x{^)d\ 



Ao 



< lim (r+i/2j^(t,Ao)) =0, 



and so (1.12) follows by (1.9). The convergence is locally uniform, since the extension 
of F\{x) I (Aa/'?/''(A2)) is jointly continuous in A > and x > (Proposition 5.4), and 
f^+^f^j^l^f^ Ao(x)) converges to locally uniformly in x > 0. □ 



Remark 6.6. Alternatively, one can deduce formula (1.12) (without uniformity in x) 
as follows. In [17j it was proved that y/tP^T^ > t) converges to V{x)l^ as t — )■ oo. 
By Theorem |1.7K a), P(r^ > t) is ultimately completely monotone. As it was observed 
in [14J, Remark 4, this already implies formula (1.12) for all > 0, see ^ (we omit the 
details). □ 



Proof of Theorem 1.8(h). The argument is similar to the proof of part (a) of the theorem. 
Again we use the notation of the proof of Lemma 6.2 Let /n,t,x(A) be the integrand 
in (Ol). Fix t > to. We have 



fn,t,x{'>^)d\ 

vWM^a(x) 



\^lj\X')dX. 



Xy/WW) 

By (|5]7|), lim^^o+ y/iAlJ^ Fx{x) / {X^/WW)) = l/r(l + a). We will use dominated 
convergence for the integral over an initial interval {0,B), and a simple uniform bound 
on the remaining interval [5, oo). 

Let B > 0. Consider x small enough, so that Ax < {it/2 — 6)/2 and iIj{1/x'^) > 2iIj{X'^) 
for A e {0,B) (recall that is unbounded). By Lemma 5.2, for A G {0,B), 



<c(e) 



'A2x2^(l/X2)^;,(l/X2) 



x,/WW) 

= c(e) 

<c(e) 

Hence, by dominated convergence. 



V''(A2)^a(A2 



'^(1/X2)(l-A2x2 
^(1/X2) -^(A2) 



< y2c(0). 



lim 



B 



1 

r(l + a) 



1/^2) / fn,U>^)dX 

) 

B 



2xNn-l/2^-tV(A2)x ,//x2 



X^'{X')dX. 
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More precisely, we have 
lim / 

x^0+ Jo 



(^/'(A2))"-i/2e-*V'{A2)A^'(A2 



T(l + a) 



dX = 0, 



and the convergence is uniform in t > to, due to monotonicity of of the integrand in 
t > to- On the other hand, 

2 



TC 



fn,t,xWd\ 



B 



which converges to as B — t- oo, uniformly in x and t > (by Remark 6.3). Hence, by 
a substitution z = t^/'(A^), 

lim (VWI^ r fn,tA>^)d\ 



r(i + «) Jo 
1 



r(n + 1/2) 



n-l/2 , = 

2r(l + a)t«+i/2 2r(l + a)t-+V2- 



Formula (1.13) follows now from (1.9). 



□ 



7. Examples 



7.1. Symmetric stable processes. These processes, corresponding to \l/(^) = |^|" with 
a G (0, 2], have already been studied in Example 6.1 in [28j. In this case i)x = (2 — Q;)7r/8, 
and Theorem 11.61 reads 



-1) — P. 



(r^ >t) = - 
with (see Example 6.1 in |28) ) 
F{\x) = sin(Ax + (2 - 0)7^/8' 



'la 



TT 



A 



na— 1 , 



e'^^"F{\x)d\ 



t, X > 0, 



(7.1) 



X exp 



la sin(Q;7r/2) 
1 



1 + - 2s°cos(a7r/2) 



TT 



1 



2a2 



log 



dc 



(7.2) 



(is. 



,0 l + C l-s°C" 
Note that all above integrands are highly regular functions (for example, log((l — 
s^C^)/(l — s°C°^)) is a complete Bernstein function), and thus ( 7.2[ ) is suitable for nu- 
merical integration. With a little effort, explicit upper bounds for numerical errors can 



also be computed. Together with (7.1 ), this gives faithful numerical bounds for P{tx > t) 
and its derivatives in t. Plots of cumulative distribution function and density function of 
Tx obtained using this method are given in Figure [Tj 

7.2. Processes with power-type Levy-Khintchine exponent. In this example, we 
assume that the Levy-Khintchine exponent ^1/(0 has the form \E'(^) = ip{C,'^) for a com- 
plete Bernstein function ip, and \& is regularly varying of positive order > at zero, 
and of positive order aoo > at ±00 (see the first part of Table [T] for some exam- 
ples). Clearly, in this case is regularly varying of orders ao/2 and aoo/2 at and 
00, respectively. Hence, by Karamata's theory of regularly varying functions (see [8]), 
iO / i^' iO converges to 1 — ao/2 and 1 — aoo/ 2 as ^ — )■ 0"*' and ^ — )■ 00, respectively. 
For simplicity, we assume in addition that the supremum in (1.7) is less than one. This 



condition is satisfied by all examples given in the first part of Table [TJ Note, however, 
that the condition given in Remark |6.5[b) is automatically satisfied, so that our extra 
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restrictions 


^2 


2 


2 


RvowniflTi TTiotion 




IS 1 


o; 




Qi-stable 


a G fO 21 


''i IS 1 ~ 


(y 


B 


(sum of two stables) 


a /3 G fO 21 a < /? 


2 _|_ 2/a\a/2 _ 


2 


(X 


vplativistir n^-stablp 


a G fO 2) m > 


((^2 ^ - l)/3/2 


/? 


a 


(subordinate relat. stable) 


a,/3 G (0,21, a< B 




2 





variance gamma 




log(l + |e|") 


a 





geometric stable 


a G (0, 2] 


i/iog(i + i/ier) 





a 


(not named) 


a G (0, 2] 


log(l + l/log(l + l/e2)) 








(not named) 




e/{i+e) 


2 





(compound Poisson) 





Table 1. Some Levy-Khintchine exponents ^(^), regularly varying both at 
(of order oq) and at ibcxD (of order aoo)- Names of corresponding subordinate 
Brownian motions are given in column Xt- First part of the table contains power- 
type functions more singular examples are given in the other parts. 



assumption can be easily dropped by referring to an improved version of Theorem 1.7 
alluded to in Remark 6.5[ b). 

Theorem |1.7K b) yields the estimate 

< (-1)" ^ P(r. >t)< ^^ i"^ 

t«+l/2^^(l/a;2) ' dt^ ^ ' - t«+l/2^^(l/x2) 



for n > 0, X > and t > Cs{n,a)/ilj{l/x ). Furthermore, by Theorem 1.8 



lim VV^(l/x2)— P(r, >t) 



-l)"F(n + l/2) 1 



7rF(l + «oo/2) t"+V2^ 
-l)"F(n + 1/2) 



TT 



V{x) 



t > 0, n > 0, 
X > 0, n > 0. 



7.3. Slowly varying Levy-Khintchine exponents. When the Levy-Khintchine ex- 
ponent "^{C) has the form '^{^) = ipiC,"^) for a complete Bernstein function ip, and \1/ is 
regularly varying of order at zero, and of order at ±oo, but at least one of Oq, Ooo 
estimates of the distribution of become more delicate. In this case Theorem 1 1.7 



IS zero, 



cannot be applied, and one needs to refer to either Corollary 6.2 (which is fairly straight- 
forward, but typically yields sub-optimal results) or the technical Lemma 6.2 Note that 
only n > 1 need to be considered, as n = was studied in general in [29 j . 

Suppose that |V'"(OI/V''(0 < 1/^ for all ^ > (that is, ^ = 1 in Theorem \L7[h)); 
this condition is satisfied by all processes in the second part of Table [T| Then = 
suP;^~^q'(9a < 7r/4 by Proposition 4.3, Furthermore, by integration, 'ip' {^i) / 'ip' {C,2) < ^2/^1 

"lib)) 



Furthermore, by integration, i^' {^1) / ^' {^2] 
when < ^1 < ^2? and therefore (cf. the proof of Theorem 



/n(t,Ao) < 



TT 



(7.3) 
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where r{a;z) is the upper incomplete gamma function. However, X^ip^ X'^) is no longer 
comparable with iP{Xq). Nevertheless, we can combine (7.3) with Lemma 6.2, to find that 



Ci(n,7r/4) 



t"+l/2^^(l/a;2 



< 



-ir^p(r.>t) 



< Ci(n,7r/4) + 



C2(n,7r/4) 



1 



t"+l/2^^(l/a;2) 



provided that t'ip{l/x'^) > 1 (so that the minimum in (6.4) is l/(t"+^/^('?/'(l/x^))^/^)) and 

(7.4) 



^Ij{1/x^ 



Here Ci{n,Q), C2{n,Q) are the constants of Lemma 6.2, and Ao(a;) = tt/{8x). 

Note that T{n + 1/2; z) < c(n)e~^/^ for 2; > for some c(n). Furthermore, ?/'' is 
a decreasing function. After some simplification (we omit the details), this gives the 
following sufficient condition for (7.4): 

ttPil/x'') > c'in) (1 + log (1 + ^(^/^') 



V V [1/x^)iIj'{1/x^ 

for some c'{n). For convenience, we state this as a separate result. 

Proposition 7.1. // |^"(OI/V^'(0 — /^^ ^ > 0, then there are positive constants 
Ci{n), C2{n), C3(n) such that 



ci{n) 



t"+l/2^^(l/a; 

/or n > 0, X > and 



<(-l)"^P(r.>t)< ^42L= 

2^ - ^ ^ dt^ ' - t«+i/2^^(i/a;2 



t > 



V'(l/a;2 



1 + log 1 



V'(l/a;2 



;i/x2)7/>'(l/x2) 



□ 



(7.5) 



(7.6) 



For example, when '^{C) = log(l + \^\") (a G (0,2]), which corresponds to geometric 
stable processes, condition (7.6) reads 

t > ^fe^ (1 + log(l + log(l + 1/x))). 



V^(l/x2 



In this case by Theorem 1.8 we also have 



lim 



log(l/x)— P(r.>t) 



Jn 



-i)"r(n + i/2) 1 



dt' 



-l)"r(n + 1/2) 



TT 



\/(x) 



t > 0, n > 0, 
X > 0, n > 0. 



7.4. Exponents with very slow growth. Let ipiC) = log(l + 'C") be the complete 
Bernstein function studied in the previous example, and let ip = ip o ... o ip be the A^- 
fold composition oi ip {a & (0,2], N > 2). The function ip is the Laplace exponent of 
the iterated geometric a-stable subordinator, and \E'(0 = is the Levy-Khintchine 

exponent of the corresponding subordinate Brownian motion. It is easy to verify that ( 1.7 ) 
is satisfied. However, (1.8) holds only if = 2 and t > 1/a {n arbitrary). Hence, for 
N = 2 asymptotic expansion ( 1.13| ) in Theorem |1.8[ b) is valid only for t > 1/a (part (a) 
of the theorem holds for all x > 0). When A^ > 3, neither part of Theorem 1.8 applies. 

Note that for any N >2 and a G (0,2], supg>o(^lV'"('C)l/'^'(0) > these examples 
do not fit into the framework of Theorem |1.7[ b), or even that of Subsection 7.3 Some 
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estimates of the derivatives of P(r^ > t) for N = 2 and t > 1/a can be obtained using 



Lemma [6. 2 1 see e.g. Corollary 6.4 



7.5. Compound Poisson process with Laplace distributed jumps. Let '^{^) = 
^^/(l + ^^). Then the corresponding process Xf is the compound Poisson process with 
Laplace distributed (i.e. with density function (l/2)e~'^') jumps occuring at unit rate. It 
was proved in [29j that 



P(tx > t) ^ min ( 1, \l —j— | > x > 0, t > 0. 



Note that (L7) is not satisfied, so the main results of the present article cannot be used. 
Indeed, we have i^x = arctanA (see \^E], Example 6.6), and so sup;^>Q is indeed equal 
to 7r/2. 



7.6. Irregular example. Estimates for are critical for Theorems L7 and L8 We 
have i^x G [0,vr/2), and it is easy to construct examples for which limx^oo^x is any 
number in [0,7r/2]. On the other hand, it seems unlikely that liminfA_s.o+ "^^a > 7r/4 is 
possible. Below we show a rather irregular example for which limsupA^o+ ''^a = 7r/2 (but 
liminfA^o+ ^a = 0). 

Consider the complete Bernstein function 

where, for example, pk = l/k\ and = 1/(A;!)^. Fix g > 0. We consider A = (qanY^'^ 
and let n — i- oo. We have 

Pk (1 + g)a„ , ^ Pk + q)an 



n-1 



< 



Pk_ 1,1 -r q)u.n ^ sr^ 
frtPn au ^£^^p„ qan 

n-1 J . oo , 

= 1 + 9 ^ + — - Yl u 

k=l ^ k=n+l 

< 1 + 9 1 1 + 9 



n 

k: 



-(n-t-l)! g(„ + 1) ^/-.^ (fc - „ - 1)! 

< (i±^ = Odin). 

qn 

Here and below the constant in the O(-) notation may depend on q. In a similar manner, 

Pk {l + qfanttk 



^^^Pn {ak + qanY 
^ Y^Pk (1 + 9)^Qn ^ Pfc (1 + g)^ 

<^^^0(1M 
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and 



o<iv/v.)i(p.(^)"-i^E 



Pk (1 + qfajak 



Pfc (1 + qfan , Pfc (1 + g)' 



It follows that 



i^{qan)=(pnY^^ (l + 0(l/n)), 

|V^"(gan)| = {pn (TT^) (1 + 0{lln)). 
With the notation of the proof of Proposition 4.5, for = ga„ we obtain 

a2 = - 1 = g(l + 0(l/n)). 



By Proposition 4.5 



limsup'i9(g„^)i/2 < lim arctan \J q{l + 0{l/n)) = arctan ^/g, 

n— >oo n— s>oo 

and, in a similar manner, 

liminft?/ -|i/2 > — I ( arcsin 




. 1-g 

= arcsm = 

4 TT 2 1 + g 

Hence, lim„_j.oo '(9(g^^-)i/2 = arctan g^/^, and therefore any number in [0,7r/2] is a partial 
limit of -(9^ as A — 7- 0+. 
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